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ABSTRACT 


A sate Hite -to -aircraft communication link over the ocean may 
include a specularly reflected signal and a diffusely scattered signal, 
as well as the direct path transmission. This report provides a computer 
program for estimating two statistical properties of the diffuse signal, 
namely the Doppler spectrum and the delay spectrum. A convenient 
representation of the total signal is also provided. Examples of the 
Doppler and delay spectra for a number of representative link geometries 
are given to illustrate the features of the diffuse signal. 
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I. INTRODUCTION 

One of the design liini tations[l ] in the development of communication 
links between satellites and aircraft is the presence of a secondary 
signal reflected from the surface of the earth. It is the purpose of this 
report to provide a convenient method for estimating the statistical 
properties of this multipath signal [2,3], particularly the power spectrum 
and delay spectrum, for ocean surfaces. Because the scattering properties 
of the sea are rather complex, it is not possible to obtain closed form 
expressions for these statistics. Instead a computer program, and its 
output for a number of representative link geometries, have been provided. 
These results should provide an insight into the characteristics of the 
multipath signal. In addition, they may be used to construct a represen- 
tation of the total signal which should be adequate for most detector 
design studies. 

The heart of the report is the computer program in Appendix I. 
Although conventional designations for range, angle, frequency etc. are 
used in the body of the report, a number of the more important identifi- 
cations between report variables and the corresponding program names are 
provided. The latter are given in capital letters enclosed in parentheses; 
for example the distance R dir from aircraft to satellite is called (REF) 
in the program. 


II. THE MULTIPATH MODEL 

Before giving the details for estimating the statistical properties 
of the multipath signal, it will be useful to describe the overall geometry 
of the link. This will permit an understanding of the mechanisms which 
produce the total receiver signal, without having to consider such 
extraneous factors as antenna patterns and polarizations. The basic 
problem in modelling the link is to estimate the relative amplitudes. 
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statistics, and signal distortions in the three components of the total 
signal transmitted from the satellite to a moving receiver. 

The first component (see Fig. 1) is the direct transmission from 
transmitter to receiver. If a voltage reference level V Q is chosen such 
that the far field radiated by the transmitter has an electric field 
strength V q /R at distance R from the transmitter (i.e., the transmitter 
power, gain etc. are subsumed in V Q ) then the receiver voltage due to the 
direct component will be, at time t 

O) »*,.(*> ■ V o< t - R di/ c > V R d1r(‘) 

where c = velocity of electromagnetic waves 

Rj.j r = TP = distance from transmitter to receiver 

and where h^ (the effective height of the receiving antenna) accounts for 
the receiver properties, including pattern and polarization, in the 
direction of the transmitter. The distortion of the signal (the Doppler 
shift for a cw signal) is accounted for by the change of R dir with time 

(2) R dl> (t) = R di> (0) + (v*n d )t 

where v is the constant aircraft velocity and n d is a unit vector along 
TP. 


The second component is that reflected from the surface (as if the 
sea were a smooth sphere) at the specular point 0. This component is 
given by 


(3) 


%ec<*> * 


V o (t - R spec /c)p 


-(khcose . )' 

1 




where 
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R = TO + OP 
spec 

= R ls + R 2s 

= distance along specular path 
p = voltage reflection coefficient at surface 
k = 2ir/X where X is electromagnetic wavelength 
h = root mean square height of ocean surface 

0.. = angle of incidence at the specular point 

h$p = antenna effective height. 

The exponential factor (the "Rayleigh roughness factor") accounts 
for the fact that since the sea surface is not smooth, some of the energy 
is scattered out of the specular signal. The distortion of the specular 
signal is also produced by the change in R S p ec with time. If the satellite 
is much further from the surface than the aircraft, 

(4 > W° ■ Vc«» + <H P >‘ 

n sp = uni t vector from 0 to P. 


Clearly the direct and specular terms are coherent, albeit with 
differing delays or doppler shifts, so their sum V^.^t) + V spec (t) will 
fluctuate (exhibit "interference", or "height gain effect") with time[3]. 


Finally, the diffuse or "multipath" component is a superposition 
of reflected fields from suitably oriented "glint points" on facets on 
the surface. If the location (which determines the delay) and the 
curvature (which determines the amplitude) of each properly oriented facet 
were known, one could construct the diffuse signal by superposition of the 
reflections from each facet. The actual sea surface is random, however, 
so that one must in practice be content with a representation of the 
diffuse signal based on certain statistical properties of the ocean surface. 
The simplest such surface is found by dividing up the mean surface into 
convenient areas A. (see Fig. 1) and associating with each area a reflected 
signal whose time delay is determined by the center point, and whose 
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amplitude is determined by the average scattering cross-section of the 
area. Thus the diffuse component can be represented by 


(5) 


diff 


(t) = I 

j 


V o (t - (R lj^2J )/c) 


1,. R~. 
lj 2 j 


a .A. 
QJ 3 
4ir 


h. 

J 


where 


R-j . = IS distance from T to scatterer 

j - SP distance from P to scatterer 

h. = receiver effective height for j-th direction 
•3 

a . = average cross- section per unit area [4] of j-th 
0J sub-area. 


For any given choice of sub-areas, each individual term in Eq. (5) will 
be in a coherent relationship with the direct and specular contributions. 
Because the range terms and R^. are continuously changing, the total 
diffuse signal will have a random character with well defined statistical 
properties (power density vs frequency, correlation coefficient, power 
density vs delay etc.). 

The sum of the three terms given by Eqs. (1), (3) and (5) now 

provide a direct representation of the total receiver signal. Insofar 

as the cross-sections a . are known, and the areas A. are appropriately 

vj 3 

chosen, the representation will possess the same statistical properties 
as the actual multipath signal. It is difficult to estimate the minimum 
number of terns in Eq. (5) needed for a statistically valid representation, 
and this problem is not considered here. In practice, the sizes for A. 

V/ 

are probably most conveniently determined empirically, by reducing their 
size and increasing their number until the statistics of interest show 
no significant change. In this report, we will not make use of this 
direct representation of the total signal, since we calculate only two 
simple statistics of the multipath signal. 
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III. THE SYSTEM GEOMETRY 


In order to introduce the program for estimating multipath signals, 
it is necessary first to establish a precise system geometry. This 
geometry is fixed by the positions of the transmitter T, the moving 
aircraft P, and the spherical earth, but is a good deal easier to 
visualize than to specify algebraically. The basic coordinate systems 
are most conveniently defined (see Fig. 2) by the specular point 0, the 
center of the earth C, and the positions P and T. The specular point is 
that point on the mean surface of the earth in the plane TPC for which 
the angle Q^(THI) between the local vertical CO and the line of sight OP 
is equal to the angle between CO and 0T. An xyz coordinate system with 
origin 0, with z-axis along the local vertical, and with x-axis in the 
POT plane may now be constructed. The tangent plane at the specular 
point is then the xy plane. The aircraft is at altitude H above the 
tangent plane with coordinates 

Xp = Htan0.j 



The satellite has range R^ s (RSAT) from 0, with coordinates 
x s = - R ls sin 6l 
y s ■ 0 

z 5 = R^coss,. 

The scatter point S on the surface of the earth has coordinates x,y,z, or 
may be defined by the polar coordinates Re,e,<f> (where Re (RE) is the 
radius of the earth) in the coordinate system x'y'z' with origin at C. 

In the computer program, variables X, Y define the scatter area center 
point at S, and the polar coordinates (STHR) = sine = (x^+y^V^/Re and 
(SPH) = <f> = tan~\v/x are computed. The distance (ZP) = (-z) = Re (1 -cose) 
of the J-th scattering point below the tangent plane, and the ranges 


6 




Fig. 2. Multipath geometry. 
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( 6 ) 


R ld = { RS ( J ) ) = ST 
R 2j = (R(J)) - SP 
R dir = (REF) = TP 

R spec = (RELSP) = T0+0P 

are now easily found. Unit vectors n p (along S to P), n T (along S to T) 
and n (the local vertical at S) have components (CNI (I , J ) ) , (CNS(I,J)) 
and (CN ( I , J ) ) along the x(I=l), y(I=2) and z(I=3) axes respectively. 

The aircraft, with position P is at the origin of a second coordinate 
system x", y", z" parallel to xyz, and is assumed to be moving with 
velocity v in the x" y" plane, with the velocity vector making an angle 
(PHV) with the x" axis. 

The antenna patterns are specified as functions of polar angles 9 , 

a 

4> a ( $ee 3) in a polar coordinate system attached to the metal of the 
aircraft, with polar angle measured from the z" axis and azimuth angle 
<f> a measured from the velocity vector. 

IV. ANTENNA PATTERNS AND POLARIZATIONS 

One must next specify the antenna patterns and polarization states 
of the transmitter and receiver. For this purpose, it will be assumed 
that T is sufficiently far from 0 and P so that the satellite antenna 
gain and polarization state is uniform over the entire scattering area, 
including the specular point and the direct path. (Otherwise it would 
be necessary to introduce still another set of coordinates centered on 
the transmitter.) Two unit vectors t T1 = - T cose, - T sine,, and T , = 

- i transverse to the line of sight may be used to specify the transmitter 
polarization via the two complex numbers A = (ATH) and B = (BPH) . More 
precisely, the transmitted electric field in the vicinity of the earth 
may be written, at range Rp 
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PROJECTION OF SP ON TO x" y" PLANE 



Fig. 3. Polarization unit vectors and antenna pattern coordinates. 



V ( t-R,/ c ) 

(7) l(x,y,z,t) = [Ai-pi+Bipi ] 

where the complex unit vector ^ 

(s) ? T = (aT t1 + bT p1 ) 

is the transmitter polarization. Since = 1» one must have AA* + 

BB* = 1. For example 

A = (0,0) B = (1,0) “horizontal" polarization 

A = (1,0) B = (0,0) "vertical" polarization 

A = (l//2,0) B = (0,+j /SZ) "right circular". 


If the transmitter is radiating a cw signal, of angular frequency «> o , or 
frequency f (F), then the reference level voltage is 


(9) 


/P t G,Z /Ait \ /2 


i 


nr ■ c- 




where Py = transmitted power 

Gy = transmitter antenna gain 
Z Q = 120tt ohms. 


The receiver antenna is specified by its maximum receiving aperture 
A , its pattern function f(e ,4> ) and its polarization state 

m cl 8 


(10) PrOVV = C<e a .* a )T T A + D (» a .* a )lp A 


where again CC* + DD* * 1 , and the unit vectors Ty^ (unit vector) and 
Tp n (unit vector) are shown in Fig. 3. In the x"y"z" system, the polar 


coordinates defining the line of sight, P to S, are e ,4> with unit vectors 

. P P 

TP’ 1 PP* 


"I TH 1 1 nn * These two sets are related by 


e a = 9 P > = > 

*a ' ^P'^V Va " 1 PP * 
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The receiver complex polarization states C = (CTH) and D = (DPH) are 
specified in terms of the polarization that the antenna would produce if 
it were operated as a transmitter. That is, the receiver antenna, if 
excited by a cw signal, would produce an electric field E^ e ^[Pp]exp[j<*> o (t-R/c)] 
in the vicinity of S where E^ e ^ is a reference level. 

Thus, for example, in direct transmission from satellite to receiver, 
the receiver voltage would be proportional to 

where 6 ',<(>' = tt-* are the values of 6 , <j> in the direction P to T of the 
a* a v a a 

transmitter, and K is a proportionality constant related to the receiver 
gain . 


The factor ^ T *^ D = AC + BD is the polarization mismatch factor. 

I K 

For example, for the receiver, 

C = (0,0) D = (1 ,0) "horizontal" 

C = (1,0) D = (0,0) "vertical" 

C = (0,j/^2) D = (l//2,0)"right circular" 

Thus the polarization mismatch factor between a right circular 
transmitter (A = 1//2", B = j/i/2) and a right circular receiver 
(C = D = 1//2) is AC + BD = j; that is, the "right circular" receiver 

completely absorbs "right circular" incident radiation. 


V. THE SCATTERING PROBLEM 

Since the intensity of the multipath signal depends on a Q , the 
scattering cross-section per unit area, we next consider scattering model 
which will provide this parameter. The predominant contribution to the 
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diffuse energy comes from a series of specular (forward) scatter points; 
the cross-section for scattering of this type at microwave frequencies 
is derived from the "optical" or "physical optics" model[4]. In this 
model, the total cross-section is the product of two factors, the first 
a polarization transforming reflection coefficient, and the second a 
probability density function for surface slope. If we consider an element 
of the mean surface of area dS = dxdy sece , at a point where the local 
vertical is n, then the geometry for specular reflection will be as 
illustrated in Fig. 4. Those facets in dS which are suitably oriented 
for specular reflection will have a local normal n 0 which is the coplanar 
bisector of n^ and np. The angle 6 ^ = (TH) - cos"l (^’^ 0 ) between the 
local vertical and the normal to the tilted specular facet then determines, 
from the known slope distribution of the surface, the probability that 
the facet normals are properly oriented. The angle between the incident 
polarization vector and the vector perpendicular to the plane of 
incidence q = n p x n T /|n p x n-J determines the reflection coefficients 
and polarization transformation. 

Thus the computational problem is to determine for each point dS, 
the voltage response of a polarized receiver due to a plane wave of known 
polarization and direction which has been reflected by a plane surface 
of complex dielectric constant e. The subroutine (COEF) determines this 
in the following manner. Consider (see Fig. 4) a coordinate system 
€,n,c with 5 the local vertical and C, n» the plane of the mean surface 
dS. A plane wave is incident from a direction 8 . in the plane. 

A scattering direction is defined by polar angles 0 scat and ♦ The 

problem is to find the polar angles 0 n <( > n which define n Q the tilted 
surface normal required to produce specular reflection from the given 
incident wave into the chosen scattering direction. The normal clearly 
lies in the plane of, and bisects the angle TSP. The scattering angle 
(the angle of incidence in the "plane of incidence" system) is o with 
2a = /TSP. The subroutine finds e n> * n and a, given e inc , e scat , * scat 
by standard trigonometric manipulation. If the incident plane wave has 
polarization ^ = aT|~| + BT p ^ and the receiver has polarization state 
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Fig. 4. The scattering coordinate system. 
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( 11 ) 


i$ 


R = C l T TS + D 1 T TS 


where the unit vectors are as indicated in Fig. 4, then it is simple to 
show, by resolving "? R into components parallel and perpendicular to the 
plane of incidence, that the effective reflection coefficient and 
polarization mismatch factor R ^ (RCOEF) is 


(12) 

II 

4 — 

O 

Cd 

R,i (a) [Acos6^+Bsin6^ ][C^cos6 s - 

•D^sin<5 



- Rj.[Bcos6 .-Asin6 i ][D-jCOs6 s +C 1 

sin6 s ] 

where 




(13) 

R,(a) 

= (eC0Sa-E)/ (ecoSa+E ) 

(RPAR) 

Ma) 

= (c 0 Sa-E)/ (c 0 Sa+E) 

(RPER) 


E - (e-sin 2 a) 1/2 


where 6. (see Fig. 4) is the angle between i-^ and the incident "parallel" 
unit vector, and is the angle between iy<- and the reflected "parallel" 
unit vector. (A complete discussion of this simple but tedious problem 
is given in Ref. 4). Thus the complete subroutine 

C0EF(9 inc- 6 scaf»scat- e - A - B - C r D l> 

accepts calculated values of the calling parameters, and returns 
a, $.j, Q n , <j> n and R^. The only problem remaining in the use of the 
subroutine is the connection between the scatter geometry of Fig. 4, and 
the system geometry of Fig. 3. By comparing these figures it is clear 
that ijj = i.pp and Tp^ = - Tpp. In the program, this is accounted for 
by changing the sign of D, i.e., C ] = C and D ] = -D (DPHI = - DPH) . Notice 
also that a slight error (entirely negligible for satellites in stationary 
orbits) is incurred in the program by assuming that the incoming wave 
polarization vectors do not depend on the scatter point position. 
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VI. THE RECEIVED POWERS 


With the coordinates and scattering parameters defined, we are 
now in a position to compute the relative power level, the doppler shift, 
and the time delay in each of the three transmission paths. 

A. The Direct Power 


The direct power received by the aircraft from a distant satellite 
is simply 

(14) p dir “ — T 21 - f(8i, »-+„) I AC + BD j 2 

air 4ir Rj. i’ v 

dir 

C = C(0. , ir-4> v ) 

D = D(6 i , ir-* y ) 


except for a small error in the mismatch factor due to the finite range of 
the satellite. The direct signal is received at a frequency f + r , 
where is the doppler shift. 


(15) 



(DDD) 


and f is the carrier frequency. The direct path doppler is printed 
out in H , a positive doppler indicating the received frequency is 
higher than the transmitter carrier frequency. In the program, the direct 
power is computed relative to a reference power level, namely that 
received by a polarization matched receiver at the specular point. That 
is, the direct power is calculated and printed as 

(16) PDIR = (R ls /R dir ) 2 I AC + BDi 2 f 

so that the actual direct power in watts is 

(in P d1r = (r T W 4 <M pDIR ) • 
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The time delay on the direct path is taken to be zero, so that all other 
delays are relative to the direct path delay. 

B . The Coherent (Specular) Reflected Power 

The magnitude of the specular path power is less than the direct 
path mainly because of the surface roughness (Rayleigh Factor) but also 
because of the reflection coefficient and the slightly longer path length. 
In the program the specular power is computed as 

(18) PSPEC- (R, s /% ec ) 2 f(’'-9 1 .''-+ v )|R cf | 2 exp(-2(klicos6 1 ) 2 ) 

where R cf = RC0EF(e.., e,. , 0, e, A, B, C, D 1 ) 

C = C (tt- 6^ , TT -4> v ) 

D-j = -D(tt- 0. , ir-* v ) 

so that again the actual specular power in watts is 

(19) P spec = (P T G T V 4 ' R ?sH pSPEC >' 

However, the value of PSPEC printed out is the ratio 10 log (P /P .. ) 

3 spec dir 

i.e., it is the ratio of specular to direct (not reference) power in 
decibels. This is a more significant ratio than the actual specular 
power for system studies. At L-band, P is negligible except for very 

o pc L 

smooth seas, or near grazing incidence. 
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c. 


The Diffuse Power 


The program computes the diffuse power by dividing the xy plane 
into squares of side dx(DX) and dy(DY) centered at x - (J + l/2}dx, 
y = (K + 1/2 )dx (J,K integers). Each square is then projected onto 
the spherical earth surface, for which the actual area becomes dxdy/cose. 
The width of the square is taken as 


dx 



2H 

2 

cos 9. 


where JJJ is an integer of order 20 to 35. This choice of square size 
ensures a smaller integration area for surfaces of small root mean square 
slope n. The incremental power scattered by this area and absorbed by the 
receiver is then given by 


( 22 ) 


P T G T A m 

dP(J,K) = - T 

4tt r:. 
1j 


f(0 ,<}> ) a 
v a ,r a oo 


cf 1 


dxdy 


4ir R2 j cose 


where a Q0 is the scattering cross-section per unit area of a perfectly 
conducting rough surface, and accounts for the reflection coefficient 
of the actual surface, and receiver polarization mismatch. The scattering 
cross-section for the sea surface, based on the physical optics model, is 
given by[4] 


(23) 


00 


= Trsece 


n 


P (0 * ) 
' n n 


where 6^, ^ are polar angles of the normal n Q required for reflection 

(see Fig. 4). Here P n ( 0 n 'l , n ) 1S the probability density function for the 

surface normal n Q , i.e., p (e n »4> n )d^ n is the probability that the normal 

lies in the cone of solid angle df! = sine de d$ . For the sea surface 

n n n n 

the actual slope distribution is a rather complicated function of the 
wind velocity and its history[5]. For the purposes of this report it is 
sufficient to assume the rather simple isotropic form 
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tan 2 e 

(24) 

P(e„ 

cos6 n n 2/ 1 . 9 „2v 

><t> ) = j e n 2n ^ (PDEN) 


-nn 

where 

2 

n = 

mean square surface slope 


= 

2 2 
<tan e n > (average value of tan e n over surface). 


This form is quite similar to a cross-wind section of the standard 
(Cox-Munk) distribution^] and, for surfaces of moderate slope, very 
nearly satisfies the normalization conditions 


(25) 


I P <Vn> di! n - 1 

) p <V„>tan 2 8 n da n - « 2 • 


Consistent with the power normalization used for the direct power, 
the normalized power increment (DP) is computed as 


/nr V 

\CO) 


D p - f<V>,)|R cf l 


•t<.n 2 0 n /n 2 (W) ^ 


Tin 


Rr,jCose 


(i) 


so that the actual incremental power in watts is 
(27) dP = DP[P t G t A m /4 T R 2 s ]. 


(It will be noted that a small error is incurred in Eq. (26) because the 

transmitter to surface distance is assumed to be R, at all times, instead 

i s 

of the variable R^. For satellites in stationary orbit the difference 
between Eq. (26) and the correct power level is insignificant.) The doppler 
and delay for this element, as before, are given by 


(28) 


f diff 


I - 
1‘ 


diff 


fn 

V 'lj 


x 2j 


n \ / _ 

"dir"*' 
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In the actual program, the incremental powers are calculated for 
four cells (at ± x, ± y) at a time, beginning with the four cells nearest 
the specular point. These four cells determine a reference power level 
(DPI). The subsequent sets are obtained by increasing |x| with |yj con- 
stant until the power from a set is less than some predetermined fraction 
of the power found in the first set. (1/1000 in the program of Appendix 
I.) At this point y is incremented, x returned to zero and the process 
repeated until further increase in y fails to produce sufficient power 
to satisfy the above condition. At this point the iteration stops, and 
the maximum values of x(XMAX) and y(YMAX) are printed. The values XMAX 
and YMAX then define an ellipse within which all significant scattering 
occurs . 

As the power from each set is computed, the doppler shift and the 
time delay for the path is also found, and the incremental power assigned 
to a doppler or delay bin. In this program 128 doppler bins cover the 
range of 2v/c f , the maximum possible range of doppler frequencies. The 
delay bins are 2 nanoseconds wide. Thus one can determine the doppler 
spectrum, (the power per Hz vs frequency) and the delay spectrum. In 
addition, the total diffuse power (PT), the sum of the power in all bins, 
is calculated in dB below the direct power. 

(29) PT = 10 log (Total Diffuse Power/P^. ) . 

VII. THE PROGRAM PARAMETERS 

To utilize the program, listed in Appendix I; a number of parameters 
must be specified. These are: 

A. Ocean Temperature and Salinity 

TEM The temperature of the ocean in °C. 

SAL The salinity of the ocean in parts per thousand. 

(Typically SAL = 35 ppt.) 

F The carrier frequency in GHz. 
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These three parameters determine the complex dielectric constant of the 
ocean via the complex function (EPS). This subroutine calculates e 
according to the formula given by Stogryn in Reference [6]. The frequency 
F should probably be limited to the range 0.4 GHz <f<5 GHz in order for 
the ocean scatter mechanism to be valid. 


B. Ocean Slope and Roughness 

ANG is the root mean square slope of the ocean in degrees. 
RUF is the root mean square height of the ocean in meters. 


In terms of the report variables, ANG = (180/7r)tan -1 n and RUF - h. These 
two oceanographic parameters may be chosen in a number of obvious ways. 
For many types of calculation however, it is convenient to relate them to 
a single variable, the wind speed W over the ocean. Thus, for example 
one may modify a standard[5] expression for the mean square slope of the 
ocean surface to obtain 

(30) n 2 « 0.006 4n(K /K ). 

n m oc 

K oc ■ 9/ M 

g - Acceleration of Gravity = 9.81 m/sec 2 
K m = Yk = y2it/X 

X = electromagnetic wavelength 
y = empirical constant - 4 


where K qc is the wave number for the largest ocean waves developed at a 
wind velocity W and is the wave number for the smallest ocean waves 
that are effective in the electromagnetic scattering process. That is, 
expressions for the mean square slope in the oceanographic literature are 
based on the actual ocean surface, including many small waves (ripples) 
superimposed on the large scale wave structure. These small scale waves 


do not participate in the forward scattering process and thus the effective 
mean square slope is smaller than the oceanographic value. The estimate 
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K m = 4 k is an empirical result which should be valid for frequencies up 
to 5 GHz and moderate to large wind velocities. Similarly, the roughness 
may be estimated from 

h 2 = .003 /K 2 c = .003 W 4 /g 2 
2 

h = .055 W /g meters 

for wind velocities greater than a few meters/second. 

C . Aircraft Velocity 

VEE - aircraft speed in meters/sec 

PHV = angle of velocity vector from x-axis in degrees. 

Thus when PHV = 0 the aircraft is flying directly away from the satellite. 

D. Geometry of Path 

THI = angle of incidence at the specular point in degrees 
H = aircraft altitude above x,y plane in meters. 

E. Antenna Polarization and Patterns 


ATH, BPH 


CTH(e * ) 
DPH{e\**) 


ANTP(e * ) 

a a 


Two complex numbers defining the polarization 
state of the transmitting antenna. 

(A, Bin report) 

Two complex functions defining the polarization 
state of the receiving antenna 
(C, D in report) 

Antenna power pattern function f(e ) in report. 


In the program in Appendix I , tie antenna power pattern f(e ,4> ) is 

a d 

taken as isotropic, i.e., f = 1. The two functions C and D are also taken 
as constants, i.e., the polarization of the receiver is independent of 
look direction. 
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F. 


Integration Parameter 


JJJ = an integer defining increment in x. 

If JJJ is taken too small, (say less than 10), the density and delay 
spectra appear erratic because the elementary areas are too large. If 
JJJ is too large, (say greater than 40) the running time becomes prohibitive. 
Reasonable values appear to lie in the range 20 < JJJ < 30. 

G. Program Output 

The output of the program, an example of which is given in Appendix 
II, are the direct power, the specular power in dB and the total diffuse 
power in dB; the specular and direct doppler shifts; the specular delay; 
and the bin widths of the doppler and delay spectra. Note that the 
unshifted carrier frequency is always in bin 65, and that the bin which 
would contain the direct and specular returns is also indicated. The 
complex dielectric constant is printed out as EPS. 

The graphical output consists of two graphs, one of power per 
doppler bin vs doppler bin number, and the other of power per delay bin 
vs delay bin number. The widths of the doppler and the delay bins are 
indicated. The total power in the doppler spectrum should sum to the 
output diffuse power level when converted to dB. 

Note that the diffuse and specular power are stated in dB "below" 
the direct power. However a negative number indicates that the specular 
or diffuse power is smaller than the direct power. Note also that when 
the receiver is blind to the direct power (right circular transmitter to 
left circular receiver, for example) the direct power is zero, so that the 
program must be modified if the specular or diffuse power is desired for 
this case. 
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VIII. RESULTS AND CONCLUSIONS 


Although the purpose of this report is to document the program of 
Appendix I, and to provide, through Eqs. (1), (3) and (5) convenient 
realizations of the multipath signal, a number of cases have been run to 
illustrate the principal features of the diffuse multipath spectra. 

In Fig. 5, the doppler spectrum is plotted vs doppler bin for 
several values of $ , the aircraft flight direction. In each graph the 
location of the specular power is shown. Evidently, the spectra tend 
to be centered on the specular signal, indicating that the major part of 
the diffuse power is emanating from a scattering region near the specular 
point. Figure 6 indicates the effect of increasing the surface roughness 
(mean square slope) on the shape and intensity of the scattered signal. 
Notice that for aircraft velocities nearly along the line of sight, the 
spectra are skewed towards the extreme doppler bins, especially for the 
rough surfaces. Note also a tendency for the diffuse power to saturate as 
the mean square slope increases. Figure 7 illustrates the variation of 
doppler spectrum as the angle of incidence is varied, and Fig. 8 shows 
the effect of polarization on the spectrum, and on the direct and diffuse 
total power. Clearly an appropriate choice of transmitter and receiver 
polarization can significantly influence the character of the total signal. 
Figure 9 shows several typical delay spectra and illustrates the fairly 
rapid initial decay and the extended tail echo for certain geometries. 

Note that in the program, the power in delay bins beyond the hundredth 
has been accumulated in the last bin. 

It may be concluded that the program of Appendix I provides a 
convenient means for estimating a number of features of the multipath 
signals which occur when a satellite to aircraft link is maintained over 
the ocean. It should be possible to undertake the study of a variety of 
detectors using either the statistical signal properties given by the 
program, or the signal realization given in Section II. In situations 
where standard detection schemes must be used, the program may be easily 
modified to investigate the feasibility of multipath rejection by 
antenna pattern or polarization control. 
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■Fig. 5. Doppler spectrum with aircraft direction as parameter. 
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Fig. 7. Doppler spectrum with angle of incidence as parameter. The oscillation in 

curve C between bins 58 and 70 is an artifact of the program, and illustrates 
the effect of choosing JJJ too small. 
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Fig. 8. Doppler spectrum with polarization as parameter. Note that the 
diffuse power in curve B is -1.90 dB below the direct power of 
0.5, or -4.90 dB below a polarization matched receiver. 
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Fig. 9. Delay spectrum for various geometries. The delay bin width is 2 
nanoseconds. A time delay of 0 (bin 1) is coincident with the 
specular signal, i.e., 14.1413 microseconds later than the direct 
signal for this geometry. 



APPENDIX I 


1 c THIS PROGRAM GIVES THE POWER SPECTRAL DENSITY AND TOTAL 

2 C *POWER IN THE DIFFUSE MULTIPATH SIGNAL FROM SATELLITE 

3 C * Tu AIRCRAFT OVFR THE OCEAN, AN j DISTRIBUTION OF POWER 

4 r * VERSUS TICE UE LAY 

5 DIMENSION SPH<4) , RS ( 4 ) ,CN ( 3 , 4 ) , cP J l 3 , 4 ) , cfJS t 3 , 4 ) . S THU 4 I 

6 n I PENSION Cl HI (4) « C THS ( 4 ) , CTHP ( 4 J , STMS ( n )« CPUS ( 4) ,SPnS(4) 

7 DIMENSION PdEL!25R ) ,DEI.(4) 

0 nif^ENS’OH PD ( 250 ) , R 1 4 ) , PHS ( 4 ) ♦ UP ( 4 ) , Ff. ; ( 4 ) ,YY (2) 

9 CO M PLl* Lp » F. HS * RCnEF * ATH « UpH,C TH ,0|'H i iH’h i 

10 50 FjRHAT(lU)', •TEfMC I = • .F6.3.5X, • G AL ( PPT ) = ♦ ♦ F6 , 3 , 5X , * F ( GHZ ) = * , F6.3, 

11 *5X» • SEfl SLOPE <nFHI = * .F5.2.5X. '.SEA HE I GHT = • , F5 , 'i « ♦ « ' ) 

12 5 1 FORMAT! 1 OX, ♦ THFTA I flt:= • , F6 . 4 , 5X , * PH I ( PL ANE VFL)='»F8.4, 

13 *5X * • PLANE HEIC-H !(«) = • ,Fl?.fe,5X, • VELOC I TY= ' , Ffl , 4 , * M/S' ) 

14 5? FOPHATOOX, ’FP<; = ',2F1().2, * O') 

lb S3 FORMAT <//.3CX« 'GRAPH OF POWER SPECTRUM B I N WIDTH I :S » , F B , 4 , • HZ • ) 

16 54 roRMAT (22X, 13.1 l*,* 12.6) 

17 b5 F0PRAr<//,10X, •01FFUSE POWER IS' .Fti.4 . •dhfcomfrent POWER IS* 

18 *,FlO, - .»* L H PF LOW OIRECT POWER OF »,F9.5.) 

19 55 FORMAT!/. ll'X.'wAX VALUES OF X AND Y IN INTEGRATION * , 5X • * XPiAXs' « 

20 *F1?.6. K X. 'Yr'AX=',El2.6. ' DElTX = '.E10.4) 

21 07 FORMAT ( ltiX . ♦ AT W= ' « 2FB .4 , 5X . «RPH= • ,2F9.4,5X, »cTH=* .2F6.4.5X, 

22 w'UPHo'.PFO.A ) 

23 56 FORMAT!//, 1(1X, .FINITE RANGE C ASE. « , 3X , ' J Jj IS *,I3) 

24 6 1 FORMAT!//, 30X, 'POWER VS DEL AY I SPECULAR r r L AY IS* ,F8, 4, 'MICROSEC ' I 

25 6? FORMAT </»30X» 'SPEC POINT IN BIN ?;bIN WXoTi! IS ',13.* WflNOSEC*) 

26 64 FORMAT!/, JOX. 'DIRECT DOPPLER IS '.Ft0.3,» HZ t 

27 +SPECULAR DOPPLFR IS '.F10.5,* HZ IN 6 I N ',13) 

26 r 

29 C 

30 C TIJPUT UAlA; TEN IS SEA TEMPERATURE IN OEc, CENTIGRADE 

31 r SAL IS SALINITY IN P.P.THOU. 35 IS TYPICAL 

32 T AND. IS RliS SLOPE OF SURFACE JN UEGRFES lg IS TYPICAL 

33 C PUF IS RMS ROUGHNESS HEIGHT PF OCEAN IN METRES 

34 c 

35 r f is carrier frequency i» gIoaherz 

36 r H IS height of aircraft above surface xn metres 

37 r VEE IS AIRCRAFT VELOCITY in methes/sec 

36 C PHV IS ANGLE OF VELOCITY FROM X AXIS)PLAnE RECEDING IF PHV=0 

39 E THt IS AfJGLF OF INCIDENCE I)' riLGRtES i 90 DEG IS grazing 

40 C JJO gives FINENESS of INTEGRATION JJOs.XO is TYPICAL 

41 f RE is RADIUS OF f ARTH IN METRES 

42 C RSA T IS UlSTANCE FROM SATFLLTTE TO SPECULAR POINT IN METRES 

43 C ATH.RPH.C1H.0PH ARE COMPLEX POLARIZATION STATESSSEE LATEP. COMMENT 

44 r 

45 c 

46 READ(6.-)TtD.SAL,F,ANG»VEE,RUF 

47 PtADlO.-JTFU.PHV.H, JJJ 

46 PEOOia.-l ATH.PPH.CTH.DPH 

49 c 

50 C 

51 P 1 =3, 1R 159265 

52 Rt = (20. /PI )*J »E6 

53 RSAT = 3.5E7 

54 WfUTeiG.bO) tem.sal.f.ang.huf 

55 WR I TE ! 6, 51 ) THT .PHV »H, VEE 
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56 WRITE <6,571 ATM,BPH»CTH,DPH 

57 wki re < 6 , bp. j jjj 

56 c 

59 C 

6,0 fta 2 =tan( Ai\iG*pr/ien. )**z 

61 C INITIALIZE STORAGE bins 

62 HO 42 IK = 1 « 250 

6 3 PU(XK)=o.O 

64 PUELUK) = 0.0 

65 48 CONTINUE 

66 P7=0.0 

67 *p=H*TAN<TH] *PT/180. I 

68 ZS = COS(THI*PT/160,)*RSAT 

69 XS = -SIiv(THI»nl/leO. )*RSAT 

70 PtF = S«RT< < XP-XS)**2 + (2S«H)**2> 

71 TDX s ? 

72 RELSP = KSAT + SQRT(Xp**? p**j) 

73 OELSP = (KELSP - REFJ/300. 

74 H2=lt*H 

75 r 

76 r. TNCIUEf'T POLARIZATION data 

77 c ATH IS THElfl COMPONENT OF SATELLITE POLARIZATION TH£ THETA UNIT 

7e c VECTOR Is -IX*rOS(THl) -l2»$TN( J t-< I ) » BPH I5 PHI COMPONENT WITH 

79 C PH I UNIT VECTOR ALONG -IVj ATH*ATHCONJ +BPH*8PHC0NJ s 1.0 

80 c VERT POL ATH=tl.P> BPh=C0,0I :HORIZ POL aTHs(O.O) BPH =<1,0> 

61 C FIGHT rlHC POL HAS ATH = (.707,0) FlPH = (0, .707) . 

82 C 

83 C CALCULATE COMPLEX DIELECTRIC CONSTANT OF SEA WATER 

84 FP = EPS(F.TF.R .SAL) 

85 WRITE (6,52) FP 

86 C 

87 C CALCUL/'Tt. INCREMENTAL AREA 

86 0TH=.S0FT(ETA2)/JJJ 

69 nX=2.*H»UTH/COS(THl*PI/iao. ) **2 

90 DY=0X 

91 XXMAXsn.il 

92 yymax=o.o 

93 c IivT E6R iTION STARTS 

94 Y=-0Y/2. 

95 I As 1 

96 3n continue 

97 Y=Y+DY 

98 IF(V.GT.YYMAX) YTMAXsY 

99 Y2=Y*Y 

100 X=BX/ii. 

101 3l continue 

102 IFIX.ST.XXHAX) XXMAX=X 

103 C TRUE POLAR COOPrmiATES OF SCATT POINT ON SPHERICAL EARTH 

104 RTHR r SLRT(X*X+Y*Y)/RF 

105 CTHR = SORT ( 1 . - STHR**2) 

106 THR = ATAN<STHR/CTHR) 

107 C ZP IS DISTANCE of SCATTER POINT BELOW TANGENT PLANE 

106 Z P = Rfc>2.*SIN( Thrt/2. )**2 

109 no 99 J = 1,4 

110 Xj = X 
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111 

112 
113 
im 

115 

116 
117 

m 

119 

120 
121 
122 

123 

124 

125 
12fe 
127 
123 

129 

130 

131 

132 

133 

134 

135 
133 

137 

138 

139 

140 

141 

142 

143 

144 

145 

14b 

147 

146 

149 

150 

151 

152 

153 

154 

155 

156 

157 
156 

159 

160 
161 
162 

163 

164 

165 


IF ( J *GL. 3) XJ = -X 
YJ = Y 

IFl J *F<>U 1 *OFU J .HQ* 3) YJ = -Y 
SPH|J) = AT AN? I Y J * X J ) 

R ( J ) = SGKTt <H + ZP>**2 + Yj**2 4* <XJ-XP>**2) 

PS < J > a SQRT< (XJ-XS>**2 + YJ**2 + (ZS+2P)**2> 

CN < 1 • J ) - SlMR*ruS«SPHCJM 
CN<2iJ> s STHP*SINCSPhC JM 
CN(3*J> = CTHR 
Ti-gKltJ) a (XP-XJ)/R<J> 

CNH2«J» a -YJ/RIJ) 

CNIC3t.J) = (H+7P)/R(J) 
fMSUi.i) a < XS-XJ)/RS( J> 

CNS<2*J) a -YJ/FSCJ) 
ri*4S43tJ) a ( 2S + 2P ) /RS ( J > 

CTHKJ) = 0. 

fTHS(J) = 0. 

C THB < J > a 0. 

SIHSIJ) = 0. 

PU 77 K a If 3 

flHl(O) = CN 4 K t J > *C NI (KtJ) + CfHlfJ) 

CTHS(J) = CNCKt J>*CNS|KtJ) + CTHS(J) 

CTHB(J) = CNIlKtJ>*CNS(Kid) + CTHB t J) 

77 continue: 

STHS(J) a SttRTll. - CTH$(o>**2) 

STHI(J) a S4RT(1* - CTHI(j)^?) 

T F t AtfSlSlHSI JW .IE* *001 .nR. ABStSTHUJI) ,LE. *001) GO TO 78 
rpHS(J) a (CTHT (J)*CTMS(J) -C T riB < J > )/(STh$[J)*STHI(J) ) 

SPHS(J) a -(YJ/APS(TJ))*StsRT<l. - CPHS<J>**2) 

PHS ( J ) = ATAN2(SPHS(J> «CPHS<JI) 

60 TO 

78 pmS { J ) =0* 

S9 CONTINUE 

f CAl-CUL'.TlON OF INCREMENTAL POWtRS STARTS 

no 40 T 1 = 1 1 4 
PHSCAT aPhSai>*160*/PI 

ThINC a ATAN<PRRT(1 # -CTHS( IT ) *#? ) /CT HS f T I ))*160./Pl 
TmSCAT = ATANCSORTCl.-CTHI ( 1 1 ) **?>/CTHl ( 1 1 ) ) *1 60 * /Pi 
C 

r THE ANTENNA PATTERN ANIP 1 1 HA ,PHA ) IS A FUNCTION Of THF POLAR COORDINATES 

C THA AN''* PPA FIXED IN THE AIRCRAFT. THE POLAR AXIS IS Z ( VERTICAL ) 

C ANO THE X** AXIS IS ALONG THF VELOCITY VFCTOP*FROM WHICH 

r. ThF anm.k pha rs measured. the ampnna POLARIZATION IS given AS TWO 
r COMPLEX FUNCTIONS OF THA.PHA. FUNCTION cTHtTHAtPHA) IS THE THETA 

r COMPONENT A^O "PH IS THF PHI COMPONENT , C TH*CTHCON J+DPH*DPHC0N J=1 

C VERT Pn|_ RtCVP HAS Ca(ltO) 0a(0*0>5 HGRl? POL RECVR HAS C a(0iO) 

C pa<liO>| RIGHT CIRC POL RLCVR HAS C=(0 t .7G7> 0=(+*707%0) 

C THA = lftU. - THlHC 

C PH A = PHSCAT - 100. -PHV 

AN TP = 1*0 

r THE NEXT LINE CONVERTS UNIT VECTORS FROh AIRCRAFT TO SCATTER COORDS 

nPHls^nPH 
C 

r ALL COEF ( TH INC • THS CAT • PHSCAT « tp • ATH t BpH f CTH * OpHl # 

* At <TH<PH<CEI I <PFLS« RCOEF ) 
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166 

167 

16 b r. 

169 

170 c 

171 

172 

173 

174 

175 r 

176 r 

177 
176 
179 
100 
181 
162 
183 c 
164 
135 
186 
167 
166 
189 r 
130 

191 

192 C 

193 

194 r 

195 
19b 

197 

198 

199 C 

200 
201 
202 

2 0 3 r 
204 r. 
20 5 r. 

206 r 

207 C 

206 r 
20 9 r 

210 C 

211 

212 r 

213 

214 r 

215 r 

216 f 

217 c 

2ia c 

219 c 

220 


nP<II) = ANTP + CA 9 S<RCOtF)**J:*PnEN<TH,PH.ETA 2 )* 0 X*UY/< 4 .*MII>*CTHR* 
*R ( II )*f 1)S(TH*PT/1P0. ) ) 

POPPLXK SHIFT IN H? PROM CARRIER IS FO* ( F*1 . E9 ) * ( V/C > 
HJ<II)^-SlN<THSCAT*Pl/len.)*COa< < PHSCAT-PHV)*PI/ 180 . ) 

DELAY IN MICRO SEC REFERRED TO OlRECT PflTrt 
ptL(II) = (R(IT)+RS< II>-RtLSP)/(300.* ( l£LSP) 

4n roNtiMtir 

TP(IA.PO.l) OPl=OP(l) 

IF(IA.EQ.l) IAsIA+1 

CALCUL * TE BIN STORAGE OF POWER* 2**7 BINS NUMBEREO 1 TO 129 
PAVG=0.« 

00 41 I J— 1 « 8 

IbS=lFlX( <1 . +pn( I J) )*64. )+l 
PUdPSIsPl ( IBS)+t)P< IJ» 

PI=PT-n>P< lj) 

PAVG=PAVb + AOS<r,P( I J> )/4 , 

CALCULATE BINS P'OR Dt.LAY 

TbT = I F i X ( 1 • + 1000,*DEL(IvJ)/lDX) + 1 

IF ( IDT .IE. 1> IRT = 1 

IF (ItlT .PL. 100) 1ST = 109 

PUCLUMT) = PDFLUPT) + DPtlJ) 

41 CONTINUE 


x=x+ox 

IF (PAVG.LT. l.E-3*0Pl. AND. X.LT. 3, *0X1 GO TO 3? 

ABOVE TEKHINATfS INTEGRATION 
IF (PAVi-.LT. 1 ,E-3*0P1 ) GO TO 30 

ABOVF. TERMINATES X INTEGRATION AND STEPS TO NEXT Y 
GO TO 51 
3p CONTINUE 

FO ri = “SIN(Thl*PT/lBO, )*C0S(PHV*P 1/160. > 

FOO = LCMPhV+PI/: f J. )*< XS-XP >/REF 
noPPLEP SHIFT IN HZ 
nut) = F{ju*F «VE F*10./3. 
cun = F(ju+f*cef*io./3 . 

IfcSrIF IX(( 1 ,+FOO) *64. ) +1 

aircraft antenna Pattern in satellite direction is 
A f 4 TP = ANTP ( Tni *PP ) 

Pp = 180. - PH V 

PECEIVFR POLAR T NATION IN SATELLITE DIRECTION IS 
CTH = CTh(THl.PP) 

OPH = OPH<THl,PP> 

ANTP=1.0 

CIRECT POWER F D ON SATELLITE 

EUlR=A'iTP*CAEiS(ATH*CTH+PPh*DPH) **?* ( PSA T /REF ) **2 

ThE RECEIVED PnWLR IN WATTS IS PnlR*POWFR*GAIN*Apf RTuRF/OIST 

WHERE PD1R HAS BEEN CALCULATED, DIST=4*PT*RSAT**2 ANO THE 

OIHER VARIABLES ARE TRANShITTERR POWE R . 1 R ANSN I T TER ANTENNA GAIN, 

PECEIVFR COLLECTING APERTURE AND RANGE FROM SATELLITE TO SPECULAR POINT 

CALCULA TL THE COHERENT POWER 

CALL COEF (THI.THI.O. ,EP, A1H,rPH,cTH,DPH1 
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221 I , f FI. $ » RCOI F ) 

222 r xl = ?.* ( <2.*Pl*F*«UF*10./3. ) *C0S(THI*PI/16C. ) >**2 

223 IF (EX1 . l~t ♦ CO.) C.X1 = fcO. 

224 FX? = EXP (-f XI) 

2?5 r>SP£C = EX 2 *ANTP*(RSAT*CA 6 S<RCCiEF)/RELSP )**2 

2?fe PS = 10.*ALOGln(FSP£C/PDIR) 

2?7 PT = 1 0 » * AL061 0 < ,J T/POIR ) 

22P. |*uOP = VtE*F/10.2 

229 PRITE<6,33) PT.PS.POIR 

230 u«ITE(b, 64)000, 000. IRS 

231 PR I TE ( h » tifc ) XXMflX, YVMAX.OX 

232 WKJTCtSt&i >V'DO° 

2 33 NN=1 

234 Y,",AX = .Ob 

235 YhIN=0.0 

2 36 no 43 I L= 1,126 

237 XX3FL0AT ( IL ) 

23b YV(l)=r>n(lL) 

239 IND=1L-1 

240 FALL PL f)T (XX ,YY,NN, IND*YMAX, YMIN) 

241 45 rO^TIi'j! ip 

242 UKlTGU.bl ) Ot'LSP 

243 UKlT£(«,b2) XOX 

244 MfJ = 1 

245 YrtAX = .05 

2 66 YfiXN = 0. 

2** 7 no 46 IN = 1 ,100 

246 XX = FLOAT UK) 

249 Y 1 ( 1 ) = POEUIM) 

250 TND = IM-1 

251 CALL P1.0T(XX,YY,NN<IND»YMAX,YMIM 

252 u 6 CONTINUE 

253 CALL EXIT 

254 FWD 
253 C 

256 c ***** SUuKCUTINES AND FUNCTION SUBPROGRAMS ***** 

257 r 

25b SUBROUTINE COFF ( THI . THS »P)iS , FP 3 » A T I , pPI , C TS t DPS , AL « Tp , PU , DELI ,OELS 

259 * • KCOE r ) 

2b0 COMPLEX ATI ,PPT , CTS, DPS .RUiEF.RPflR « RPER . EPS 

261 C GIVEN Tut. INCIDENCE AND SC A TTfc'H 1 NO DIRECT IONS, THE POl. AR 1 7 ATI ON 

262 C STATE OF THE TRANSMITTER AMP RECF I VFP , AND THE COMPLEX DIFLECTRIC 

263 C CONST A*'T t THIS SUPROUTIWE pF TrRrllNF.S THE DIRECTION OF THE NORMAL 

264 C TO THE SURFACE Rr OuIPfLi FOR ?UCH REFLECTION > THE ANGLE ALPHA 

263 c BETWEEN THF NORMAL ANO FIT HE p THF INCIDENT OR REFLECTED 

266 C DIRECTIONS, AND THE COMPLEX RFFlEcTiON COEFFICIENT, 

267 c ALL ANGLES REGiJIREO IN pE&RFFS. 

266 PI=3. 14136265 

269 r THI -THETA ANGLE OF INCIDENCE 

270 r THS-IHPTA AI&LF OF scattfr 

271 C PHS- PHI ANGLE OF SCATTER 

272 C AlI-riCTA COMPONENT OF THE INCIDENT FIELD 

273 C PpI-PHl COMFONF NT OF THE IMCIOENT FIELD 

274 c CTS-THFTA COMPONENT OF THE SCATTERED FIELD 

275 r OPS- PHI COMPONENT OF THE SCATTERED FIELD 
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276 r 

277 
276 
279 
2BU c 

261 r 

262 r 
2C» 

264 
266 
266 
26 7 
2B£ 
269 

290 r 

291 C 

292 
296 

294 

295 

296 

297 
296 
29 9 
30C 

301 

302 

303 

304 

305 
3»f> 
307 
30c 

309 

310 

311 r 

3 1 2 r 

313 r 

314 

315 
3lb 
317 
31C 

319 

320 

321 

322 

323 

324 

325 

326 

327 
326 

329 

330 


rHS-CO'"PLfeX DIELECTRIC CoNSTaN r 
ThI=THT*PJ[/lPn, 

THSs:THS*ri/J 60. 

PHSxPhS*Pj/3 60. 

AL-ANbLE ALPHA RF Tu'EtN I PC IDFNT DIRECTION A NO NORMAL 
rH- PHI ANbLf. OF the NOHWaL 
Th^THLTA AM ; LF OF THE NCRhAL 

r^AL = -CSlfi( 1 Hi >*SIN( ThS > *lOS ( P n S > -COS < T H I ) *COS ( T hS >) 

FALsSasTi ll.+CSM )/^. ) 

A L - A T A f 1 F < S U R T f l . r At) t*?) ,CAL) 

T h ( A L * * . F • A * L - 3 ' ftO TO 100 

Ph=ATAMS> l SIXTHS) +S IN<PHS) tSTN ( THS ) *COS< pHS I -SIN ( THI ) ) 

CTH=(COS( ThS)+rOS< THI l ) /(2.*CAL) 

T h * A T A o' 2 ( S 0 K‘ T ( 1 .~CTH**2 ) ,CTH) 

HtLI-A'inLL OF LT A I NClDENCt 
PtLS-A fGLL i‘li. TA SCATTERED 
■s < *AL = STNl?,*AL> 

T uFL I = < SIN < ThS)*COS<PhS ) *COS( ThI )+SIN l ThI ) *CoS(ThS > ) /S2AL 
SjFLl = t: lN (THS ) *RIN< PHS) /S<AL 
ntLIsATAN^i SOFl.lt CDf-LI) 

ruELS= (SXN( THI ) *COS( THS ) ♦COS ( F hS ) + COS C THI 1 *SI C T MS > ) /3£AL 
SuELS = SIM| THI > *SIN< PHS) /$2 AL 
ntLS = ATANif < SOELStCOELS) 
ru TO 90 
ICO CONTINUE 
A L = u • 0 
THsTHl 
PH-PI 

flliso.o 

PtLS=U.O 
CUF.LI=M 
SuELI = ').u 
CuELS^I.O 
£uELS^«i*0 
90 r u n T I n 1 1 f. 

PPAR- F'^SNtL REFLECTION COEFFICIENT FOR PA^At LEL POLARIZATION 
PF'FR-F^tbFLL REFLECTION COEFFICIENT FOR PERPENDICULAR POLARIZATION 
PC°EF-ror-f-LE.X PEFLECTlON COFFFiClE'NT 

RpARsi *- PS+COS ( aL>-CS0KT<EPS-*In( AL>**2) ) /(EPS* COS( AL) +C SORT (FPS-SI 
4 N < A L > * * 2 > > 

PPER*(rOS( At. )-rSORT<tPS-SIN( aL)** 2> >/<COSl AL > +CSQRT ( £Ps-SIN C aL > **2 
*> > 

RCOEF = -RPEK*< ( P.P I«XDEl I - ATI *SDLL I ) * ( DPS*CUEl S+C TS*SOFLS ) ) +RPAR* ( < A 
*Ti*COtl.I + '-|-I*snft I ) ♦(CTS*CnruS-DPS*SDELS) J 
A L- AL 1 .q(j . /HI 
TH=TH* lfib./Hl 
PH=PH*13U ./PI 
THI=THT*1 tfO./PT 
ThS = THS*ieC */PT 
PhS=PHS*lf G./PT 
OLLIsGELlvlA0./PI 
OLLSsDFLS*16C • /PI 
RlTURn 
Erin 

FUNCTION PDLN«THtPHtS2) 


34 



331 C This FUNCUOM CALCULATES THE VALUE OF THE PROBABILITY OENSITY 

332 C FUNCTION ( H f ! E N ) FOB A SUPFACF NORMAL WHERE TH IS THE THETA-ANGLE 

333 r IN DEiji'F.t-S ANP PH IS THF PHI-ANGLE IN DEGREES OF THE SURFACE NORMAL 

334 C AND S? IS THE MEAN SGUARE SLOPE. THE PDF IS GAUSSIAN. 

335 Pi = 3.141b«2i-5 

336 TH=TH*PI/U«'. 

33 T A = -TAN( Tl')*l AN{TH)/(S2*(1.+2.*S2) ) 

33d IF < AHS< AJ .GT.fcf). ) A=-60. 

339 PL)EN=COS<TH»*PVP( A)/(PI*S2) 

340 TH=TH*1B0./PI 

341 RETURN 

34^ F Nl! 

34 3 COMPLEX I- UNCTION EPSIF.T.S) 

34lt r COMPLhX 0 1 ELECTRIC OF Sf A WATER EPS ISfOGRYN MTTJ 

345 r T IN rrMlGRAOE. F In GIGAHFR12. s IS SALINITY in p.p.thou. 

346 kFAL n 

347 C NORMAL Sea water RANGES FROM 20. To 35. PPT 

340 N = SM. 01707 + .00001205*S + S*S*<* . QSRf -9 > 

345 C = 2*>. - I 

350 AA = 1 .6451—5 - C * ( 1 . - . l*r I «2.551E-7 

351 U = LXP(C.*< .0^033 + . Un01? b t-*^+E*r*2.464E-6-R*AAM 

35?. t,IG = S*(.l6P5?l-.Of)14615*S+.OOn02093*4*S - St S*S* 1 . 262E-7 ) /O 

353 00 =.11105 - , 003p24*T + . CO 0 0654*1*1 

354 b = ( 1 . + ,0C14F3tM*T - - . 02 e J67*N*P )*00 

355 fO = < l.-.2551*M+.0515*N*i\i)* (07.74-. 400P+T+. 0OO94*T*T ) 

356 . EPS = 4.5 4 (FO-4.9)/CMPLX( 1. » W*F> - CMPEX ( 0 • * 1 7 • 9S*S I G/F > 

357 RETURN 

350 Fi'iO 

359 SUPROUTlNtPLOTI X • Y . N . I NP . Y MAX . YM'IN) 

300 OIMENSIONrtl 1191 .YLAITEU 6 ) . Y <1 0) .MARK ( 10 ) 

3fil r>ATA MARK < 1 ) .MARK <? I .MARK ( 3) .MakK (5) . M AOK <6) .MARK (7) .MaRKIS) ♦ 

362 2MAPK ( 9 ) « Mi ARK ( 10 ) i MARK (4 ) /1H* . lh, , 1HI . 1H0. 1HN.1HH. 1H1 . 1H2. 1H- , 1HX/ 

303 OATA I !, L AMK • NOP T « If'LOS/lH .1HS.1H+/ 

364 IFdNHll. 1,11 

365 1 WRITE H> * 3 ) 

366 3 FORMAI (//25X,*ORPER IN WHICH PLOT SYMBOLS ARE USED * . IXONH1Z- • t 

367 *//30X « • THE SYMBOL 1$) INDICATES OFF-SCALE DATA**//) 

366 007J=9.H9 

369 7 K( J)=MAHK(10 ) 

370 MtOUNT=lo 

371 PCAUE = lOT..O/( YMAX-YMIN) 

372 i ll=( -YM1N*SCAI.E>+11.5 

373 hoej=i .6 

374 P=J-1 

37b 0 YLA 8 EL< J»sh*20.0/SCALE+YM1N 

376 WRITE (4.6) lYL'BtHI) .1=1*6) 

377 9 pyRMAT (6X.lPE9.2*5(lPE20,2) / ) 

37t> G0T0132 

379 11 MC0UNT=NC0L‘NT + 1 

300 P099U=1 . 119 

351 99 M<J 1=1 PLANK 

36? TF (LLL.GE.ll .AAfl.LLL.LE.no >M(lLL)=MARk (10 ) 

303 IF (NC0>iNT-10)133.132*l33 

3ii4 135 1)0890 = 11. 111. 20 

305 09 M(J)=IPUUS 


35 



3B6 

133 

ri)20J~l t<’J 

3e; 


L=< Y< J)-’rMU)*SCALr + 0,5 

3*G 


T l- { L M “ , 1 7 * 1 7 

3<vj 

1 4 

T f t L + 1 i ) 1 5 ♦ 1 6 1 1 6 

390 

15 

1 )=>H.'4 J 1 

391 


nUT02U 

392 

16 

l.lxl*u 

3^3 


MLL>-MAk|<< j) 

39* 


r.uTo^u 

395 

17 

TF (L-10fU10tl9tl9 

396 

IB 

LI = L * j 1 

397 



396 


PijTOau 

399 

19 

fM 119>r r *iUf ; T 

*4 0 0 

20 

rO^T IN' 4 : 

401 


TF <fcCU‘ [ NT-10 >?1 t2b * 21 

402 

?1 

UrtITE(--»K4> <I”(J ) f Jsl *119) 

403 

24 

r ylViA T t IX , 119A1. ) 

404 


( : j T f J 2 7 

403 

2 5 

iFhITL<^*2f- > (Xt € M C J » tj-9tll9) ) 

4"6 

26 

FU^MAT ( 1X.F6.3 1 HlAl) 

407 


NCOUNTsO 

406 

2? 

CONTINUE 

409 


RfcTUHN 

410 


F fcUS 
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APPENDIX II 

SAMPLE OUTPUT {Curve B Fig. 8) 


TtM(C)=l 0,000 PA UPPT 1*35,000 fl 6 * 17 ) = 1,400 SEA StOPE<GEG)=l5af )Q SEA HEIGHT^ ,30tf 

THP1A IMIs ^S.OOnO PhKPLAMF VfL)= 90,OOCG PLANE HEiGHT(f*iJ= , 3 U 0 DIK>E 4 VELOC IT Y = 20 U * OOO OM/S 

AlH= *?n ?0 O.nOUO RPh= O^OOOO *7070 CTH= 1*0000 0,0000 UPH= 0,0000 U,UoQO 


PIMTE CASr jjj IS 30 

F pPs; 71,57 -5f.,P3 J 


DIFFUSE POWER IS -l*tt943CBiC0HFRENT POWER IS -?62«304J OB BELOW DIRECT POWER OF .499eS 

DIRECT OORPlER -.COO H7J SPECULAR DOPPLER IS **000 HZ IN BIN 64 

WAV VALUES OF X *Nti Y IN INTEGRATION vv,AX= ,3lG2et>E b YMAXs ,bl?ti21E 4 DELTX a *1072£ 3 


GRAPH OF POWs.« SPECTRUM BIN WIDTH IS !4,5a35HZ 


o»dfr in which plot symbols ARE UStD * • I XONHIZ* 

THE SYMBOL <*) INDICATES OFF-SCALE DATA 


0 « 0 0 1 0 


11,000 * 

a 

* 


* 

* 

* 


1 * 0 0 £ .2 2,00- -2 5.00E -2 4»Q0£ -2 

+ 


+ 


*■ 


+ 


+ 


5 . 00 £ -2 

--**+--* 


+ 


21,000 


+ 


♦ 


+ 


+ 


♦ 


31,000 


+ 


♦ 


+ 


+ 
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si.noo + 


61. TOO + 


71*0(j0 + 


01 □ 0 oo + 


91,000 + * 
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101 ,000 


- ■* 
♦ + 


4 


4 


4 


♦ 


4 


111*000 


* 

* 

* 

* 

* 


121.000 * 


+ 


4 


4 


4 


♦ 


4 


4 


4 


4 


POWER VS CTtAY jSpECULAR DELAY IS 1<4 . 1 *U3MlCR0SEC 

SPEC POINT IN BIN WIDTH IS 2 NANOSEC 

OPOF.R IN WHICH PLOT SYMBOLS ARE USED *.IXONHlZ- 

th£ sywbol U) Indicates o^f^scale data 


O.OOL 0 

1,000 — »-•*-“ 


1.00E -2 

♦ 


2 . ocr -*2 


3.U0E -2 

— 


4.00E -2 

— *-- 4 ~ — - 


5« OPE -2 
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11*000 


+ 


21.000 + 


31. COO + * 


*♦ 1.000 4 * 


51. TOO + 


61.000 + * 


71,000 + * 


4 


♦ 


4 - 


♦ 


4 


4 


4 


♦ 


4 


4 


4 


4 


4 


4 


4 


4 


4 


4 


4 


4 


4 


4 


4 


4 


4 


4 
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- * 

~ * 


ei.noo 4 * 

- 4 


+ 


+ 


* 

* 


91,000 


- * 
-r * 


m * 
- * 


- * 
- * 
- * 
- * 
— * 


4 


♦ 




41 



REFERENCES 


1. Painter, J. H., Gupta, S. C., and Wilson, L. R. , "On the Technology 
of Aerospace Communication in Multipath, IEEE Nat. Telecomm. Cont., 
Houston, Texas, December 1972. 

2. Durrani, S. H. and Staras, H. , "Multipath Problems in Communications 
Between Low Altitude Space Craft and Stationary Satellites," RCA 
Review 29, 1968, p. 77. 

3. Gager, D. R. , "Signal Fading Characteristics due to Reflection From 
the Sea, Marconi Rev., First Quarter, 1971, pp. 27-48. 

4. Ruck, G. T., Barrick, D. E. , Stuart, W. D., and Krichbaum, C. K., 

Radar Cross-Section Handbook , Plenum Press, New York, 1970. 

5. Phillips, 0. M. , The Dynamics of the Upper Ocean , Cambridge University 
Press, 1966. 

6. Stogryn, A., "Equations for Calculating the Dielectric Constant of 
Saline Water," IEEE Trans. MTT, 1971, pp. 733-735. 


42 



